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We have studied the phase diagram of the one dimensional S — | XXZ model with ferromag- 
netic nearest-neighbor and antiferromagnetic next-nearest neighbor interactions. We have applied 
the quantum renormalization group (QRG) approach to get the stable fixed points and the running 
of coupling constants. The second order QRG has been implemented to get the self similar Hamil- 
tonian. This model shows a rich phase diagram which consists of different phases which possess 
the quantum spin-fluid and dimer phases in addition to the classical Neel and ferromagnetic ones. 
The border between different phases has been shown as a projection onto two different planes in the 
phase space. 

PACS numbers: 75.10.Jm, 75.10.Pq, 75.40.Cx 

I. INTRODUCTION 

There is currently much interest in quantum spin systems that exhibit frustrations. This has been simulated in 
particular by study of the magnetic properties of the cuprates which become high- T c superconductors when doped. 
Frustrated spin systems are known to have many interesting properties which are quite different from the conventional 
magnetic systems. 

The Heisenberg spin i chain with nearest neighbor (NN) and next-nearest neighbor (NNN) interactions (which is 
equivalent to a zig-zag ladder) is a typical model with frustrations. In the recent years, several interesting quasi-one- 
dimensional magnetic systems have been studied experimentally [3, 0, Q • Among them, some compounds containing 
CuO chains with edge-sharing CuOi plaquette were expected to be described by the XXZ model with next-nearest 
neighbor interactions. The nearest-neighbor (Cu — Cu) spin interaction changes from antiferromagnetic (AFM) to 
ferromagnetic (FM), as the angle 6 of the Cu — O — Cu bound approaches 90°. The next-nearest-neighbor interaction 
is always AFM and is not dependent on 6 Q. Several compounds with edge- sharing chains are known, such as 
Li 2 Cu02, La§C a%Cu2\0 C^i^C^sOio, Rb2Cu 2 Mo 3 Oi2, which can be considered as an ideal model compounds 
with the ferromagnetic NN interactions and antiferromagnetic NNN interactions Q . 

The Hamiltonian of such model on a periodic chain of N sites is 

N N 

H = 4 { E« CT ?+i + + A CT ^i) + E -M« +2 + ° y y i+ 2 + of <+a)>, (!) 

i=i t=i 

where J > and J2 > are the first and second-nearest neighbor exchange couplings and the corresponding easy- 
axis anisotropies are defined by A and A 2 = J26. For J 2 = 0, the ground state properties are well known from 
the Bethe ansatz Q- For positive coupling constants (J, J2,A, 5 > 0), this model has been investigated previously 
[H, Q. In particular, it has been shown that a transition from a gapless state to a dimerized one takes place at 
( J2 = 0.24, A = 5 = 1) 10]. The point (J2 = 4, A = 6 = 1) corresponds to the well known Majumdar-Ghosh model 
where the exact ground state is constructed from the direct products of dimers which leads to a gapful phase [Til ] . 
Relatively, less is known about the model with the ferromagnetic NN and the antiferromagnetic NNN interactins. 
Though the latter model has been a subject of many studies O 0, 0, EE] , the complete picture of the phases in this 
model is still in investigation [l6j |. It is well known that there is a critical point (J2 = 4, A = —6 = —1) where the 
ferromagnetic state is unstable and the ground state is nontrivial at J2 > \ which can be realized by different phases 
[l7j |. Moreover, the exact ground state can be represented in the resonating valence bound state (RVB) [11,0]. This 
state has been proposed as a candidate for the spin liquid ground state [201 ] . One of the most important and open 
question is the possibility of the spontaneous dimerization of the system in the singlet phase accompanying by a gap 
in the spectrum [2l|. The controversial conclusion exists about the presence of a gap at J2 > \. It has long been 
believed that the model is gapless [12, but one loop renormalization group analysis shows [3, [24| that the gap 
opens due to a Lorentz symmetry breaking perturbation. However, the gap has not been checked numerically [14j . 
On the base of field theory consideration it was proposed [25| that a very tiny but finite gap exists which can not be 
observed numerically. 

In this paper we have considered the one dimensional anisotropic S = \ Heisenberg model with ferromagnetic NN 
and antiferromagnetic NNN interactions by implementing the quantum renormalization group (QRG) method. We 
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FIG. 1: The decomposition of chain into three site blocks Hamiltonian (H B ) and inter-block Hamiltonian (H BB ). 



have calculated the effective Hamiltonian up to the second order corrections. The second order correction is necessary 
to get a self similar Hamiltonian after each step of QRG. In this approach, we have considered the effect of whole 
states of the block Hamiltonian which are partially ignored in the first order approach. The present scheme allows 
us to have the analytic RG equations, which give a better understanding of the behavior of system by running of 
coupling constants. We have succeeded in obtaining the phase diagram in a good qualitative agreement with the 
numerical ones [26|. 

We have previously studied the antiferromagnetic model (Eq.([T])) for A > by QRG |8j. For < A < 1 the 
interplay of the two competing terms (NN and NNN) in the presence of quantum fluctuations produces the dimcr 
phase for Ji > J|(A, A2). The dimer or spin Peierls phase has a spin gap and a broken translation symmetry (the 
unit cell is doubled) in the thermodynamic limit. However, we have determined the fluid-dimer phase transition by 
using the running of couplings under RG (see Fig. 3 of Ref.@ or the complete phase diagram presented in Fig [3] in 
this article). In the spin- fluid phase, the anisotropy and next-nearest neighbor couplings are irrelevant while in the 
dimer phase they run to the triple point (Aj = J| — 0.155, A* = 1). From a quantitative point of view at A = the 
RG analysis gives J| — 0.44 which can be compared with the numerical result of Jf ~ 0.33 presented in Ref.[9(. The 
Neel phase appears just by crossing the A = 1 plane at A2 = and J2 = 0. In the A2 = plane and for A > 1, the 
model will pass through a phase transition from Neel to dimer phase for J2 > J|(A). The Neel ordered is also broken 
by increasing the anisotropy of the NNN interaction. 

In this paper we will complete the phase diagram of this model for the whole range of parameters. Moreover, we 
intend to consider the XX Z model with ferromagnetic NN (J < 0) and antiferromagnetic NNN ( J J2 > 0) interactions 
which can be fulfilled by extending the phase diagram to J < and JJ2 > 0. If we implement a 7r rotation around z 
axis for the even sites and leave the odd sites unchanged, the Hamiltonian (with J < and JJ2 > 0) is transformed 
to the following from 

j N N 

H = I { E«<+i + ° V i°U ~ A °f°?+i) + E -M« +2 + + 5 ^< +2 )}; J,J 2 ,A,S>0. (2) 

i=l i—l 

The QRG procedure is implemented on the rotated Hamiltonian (Eq.([2])) which makes the calculations easier. How- 
ever, the phase diagram and other figures presented in this article are based on the couplings defined in Eq.([T]). 

The QRG approach will be explained in the next section where the second order effective Hamiltonian and the 
renormalization of the coupling constants are obtained. In Sec. Ill , We will present the phase diagram and its 
characteristics where a comparison with numerical results is done [26| . Finally, we summarize our results. 

II. RG EQUATIONS 

The main idea of the QRG method is the mode elimination or thinning of the degrees of freedom followed by 
an iteration which reduces the number of variables step by step until a more manageable situation is reached. We 
have implemented the Kadanoff's block method to do this purpose, because it is well suited to perform analytical 
calculation in the lattice models and they are conceptually easy to be extended to the higher dimensions. In Kadanoff's 
method, the lattice is divided into blocks where the Hamiltonian is exactly diagonalized. By selecting a number of 
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low-lying eigenstates of the blocks the full Hamiltonian is projected into these eigenstates which gives the effective 
(renormalized) Hamiltonian. The effective Hamiltonian up to second order corrections is [U, [U HH 

Jjeff = jjeff + H eff + R eff 

H* ff =P Q H B P a , H{" =P H BB P Q , Hl Sf = P [H BB (1 -P$ 1 nB (1 - P )H BB }P . 

We have applied the mentioned scheme to the Hamiltonian defined in Eq.([2]). We have considered three-site block 
procedure defined in Fig.([T]). The block Hamiltonian (H B =^2hf) of the three sites, its eigenstates and eigenvalues 
are given in appendix A. The three site block Hamiltonian has four doubly degenerate eigenvalues (see appendix A). 
Po is the projection operator to the ground state subspace which defines the RG procedure. Due to the level crossing 
which occurs for the eigenstates of the block Hamiltonian, the projection operator (Po) can be different depending on 
the coupling constants. Therefore, we must specify the regions with the corresponding ground states. The eigenvalues 
of the block Hamiltonian are labeled by eo,ei,e2,e3 (see appendix A). In the following, we will classify the regions 
where each of this states represent the ground state. A summary of this information is given in Fig. (J4J) of appendix 
A. 



A. Region (A): eo is the ground state. 

In this region the effective Hamiltonian in the first order correction leads to the XXZ chain without the NNN 
interaction (J 2 — 0), i.e the effective Hamiltonian is not exactly similar to the initial one. The NNN interaction is 
the result of the second order correction. When the second order correction is added to the effective Hamiltonian, 
the renormalized Hamiltonian, apart from an additive constant, is similar to Eq.([5]) with the renormalized couplings. 
Thus, the effective Hamiltonian including the second order correction for A > is: 



H eff = _ 



J 1 



N/3 N/3 



The renormalized coupling constants are functions of the original ones which are given in appendix B. 

B. Region (B): ei is the ground state. 

The second order effective Hamiltonian is similar to the case of region A with different coupling constants given in 
appendix C. A note is in order here, although the second order correction is necessary to produce the NNN interaction 
in the effective Hamiltonian the initial values of J 2 = and A2 = do not produce NNN interactions. It is different 
from the RG flow obtained in region A. 

C. Region (C): e$ is the ground state. 



In this region the effective Hamiltonian to the second order corrections leads to the Ising model 
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FIG. 2: The cross section of the three dimensional phase diagram for A < 0. (a) The crossed plane is A2 = J2A and (b) 
A 2 = 0. 



This simply introduces the ferromagnetic behavior. We will discuss the phase diagram in terms of different regions 
defined above in the following section. 



III. PHASE DIAGRAM 



A. Region (A) 



In the A < case, the RG equation shows running of J to zero which represents the renormalization of the energy 
scale. We have plotted the RG flow and different phases in Fig.©. The solid line is the boundary between dimer 

(II) phase and spin-fluid (I) phase. If we start from the spin -fluid (I) phase or Dimer (II) phase, the sign of NN 
anisotropy (A) changes under RG after few steps. However, in the dimer (II) phase the amount of NNN coupling 
( J2) is greater than 0.44 just when the anisotropy changes sign, while in the spin-fluid (I) phase it is less than 0.44. 
In other words, in dimer (II) phase the RG flow goes to the triple point (A| = J| — 0.155, A* = 1) (the filled circle 
in Fig. ([3])) while it goes to A = J2 = A2 = fixed point starting from the spin-fluid (I) phase. In Fig. ©a)) and 
Fig- (Oh)) the black arrows show the running of couplings under RG. In the region denoted by (g), though eo is the 
ground state, the behavior of the couplings constant is not the same as the couplings in the dimer (II) phase. In this 
region ((g)) coupling constants go to the spin-fluid (III) phase. It means that the region denoted by ® and spin-fluid 

(III) are a unique phase. The authors in the reference [261 ] were not able to specify the phase of this region numerically. 
We denote the boundary between the spin-fluid (III) and both dimer (II) and spin-fluid (I) phases by long-dashed 
line. The dashed line behind the arrows is not a phase boundary and just represents the two regions with differenet 
ground states (see appendix A). It is known that on the A2 = J2A plane and A = —1 there is a fixed point, namely: 
A = -1, J 2 = 0.25 DHHi. 

However, our approach is not able to show this fixed point because this is on the plane 
which is separated by spin-fluid (I) and ferromagnetic phases where the level crossing occurs. Instead, close to the 
A = — 1 line we found the critical value of J| = 0.28 which distinguishes the spin-fluid (I) and spin-fluid (III) phases. 



B. Region (B) 



In this region the NNN interactions are greater than NN interactions. The implementation of bosonization technique 
combined with a meanfield analysis in the reference [24[ predicted that for A2 = J2A plane and A = 0, the system 
might exhibit a chiral ordered phase with gapless excitations where J2 > </f = 1.26. The predicted critical value (JJ) 
is in well agreement with the numerical density matrix renormalization group result (30| . Our approach shows that all 
coupling constants are irrelevant except J2 and A2. For — 1 < A < the ratio of A2 to J2 goes to zero and for A < — 1 
this ratio goes to infinity. It means that in the fixed point of this region the original spin chain decouples to two XX Z 
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FIG. 3: The projection of the complete phase diagram on the A2 = J2A plane in (a) and on the A2 = plane in (b). 



chains without next-nearest-neighbor interactions where the lattice spacing is doubled. For — 1 < A < 0, the model is 
in the spin-fluid (III) phase which is specified in Fig. ([2]). The spin-fluid (III) is different from the spin-fluid (I) phase 
according to their stable fixed points. The stable fixed point for spin-fluid (I) is J 2 = 0, A2 = while for spin-fluid 
(III) it is ^ — > having J 2 , A 2 — > 00. Note that the level crossing of eo and e 2 does not define the border between 
spin-fluid (III) and dimer (II) phase. This border is defined by the running of couplings under RG equations. For 
A < — 1 the model is in the antifcrromagnetic phase. In this case the model is decoupled to two antifcrromagnetic 
Ising chains. Thus the ground state is long-ranged antiferromagnetic ordered ITT44TT44''')- 



C. Region (C) 

As we pointed out in sec[Tl]-C, even after adding the second order corrections, the original Hamiltonian is mapped 
to the ferromagnetic Ising model. Ising model remains unchanged under RG as fixed point and its properties are well 
known. We call this region as the ferromagnetic phase. 



IV. SUMMARY AND DISCUSSIONS 



We mapped the one dimensional ferromagnetic NN and antifcrromagnetic NNN S = i XXZ model to the anti- 
ferromagnetic model in Eq.((2]) with negative anisotropy. We have implemented the second order QRG procedure to 
get the phase diagram of this model. The complete phase diagram which also covers the positive anisotropy region is 
presented in Fig. ([3]). This is a cross section of the phase diagram with A 2 = J 2 A plane in Fig[3ja) and with A 2 = 
plane in FiglHJb). For —0.5 < A < 1 (on A 2 = J 2 A plane), when J 2 is smaller than the critical value (Jf) the system 
is in the gapless spin-fluid (I) phase. By contrast, for larger value of J 2 > J 2 the system is in the dimer phase with a 
finite energy gap above the doubly degenerate ground states (the transition is denoted by solid line). As J 2 increases, 
the system exhibits a transition from the dimer phase to the spin-fluid (III) phase which is called the gapless chiral 
phase in Ref.[30(. The transition is denoted by long dashed line on the phase diagram. The transition takes place 
at J| = 1 for A = 0, in qualitative agreement with Jf = 1.26 of Ref.^. For -1 < A < -0.5 (On the A 2 = J 2 A 
plane) at J| a transition occurs from the spin-fluid(I) to spin- fluid (III) (chiral order) phases. The QRG equations 
for A > — 1 shows a critical line (long dashed line) which separates the spin-fluid (I) and spin-fluid (III) phases 
without an intermediate region. Thus, we claim that for J 2 > Jf = 0.28 there is no gap and the model is not in the 
dimer phase. The model is in ferromagnetic phase where A < — 1 and small J 2 . The phase transition to long-range 
antiferremagnetic phase takes place at J 2 > J% (dashed-dot-dot line). In the case of A > 1, a transition from Neel 
(VI) phase to the dimer phase occurs as J 2 increases (solid line) . The dimer phase is unstable by increasing J 2 further 
which leads to a transition to the antiferromagnetic (AF(III)) phase (long dashed line). However, The comparison 
of FigGJa) with FigJ3]Jb) shows that the anisotropy of the NNN-term (A 2 ) changes the phase diagram significantly. 
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From the parameters estimated for several compounds near the isotropic limit (A = 1, A2 = J2A), La§Ca%Cu 2 \0^\ 
(Jf = 0.36), Li 2 Cu0 2 {J 2 = 0.62), Ca 2 Y 2 Cu 5 O w (Jf = 2.2) Q, our result predict for A > -1 that all of them are in 
the chiral order phase (spin-fluid(III)) without gap, and for A < —1, La§Ca%Cu 2 \0±\ is in the ferromagnetic phase 
and Li 2 Cu0 2l Ca 2 Y 2 Cu^O\Q are in long-range antiferromagnetic order. 
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APPENDIX A: THE BLOCK HAMILTONIAN OF THREE SITES, ITS EIGENVECTORS AND 

EIGENVALUES 

We have considered the three-site block (Fig.fl])) with the following Hamiltonian 

hf = 1 Mjvli + oljolj + al jal j + a v 2J a v 3J ) - A(oljO z 2J + o^alj) 

+H°i,i°li + <i°h) + A 2 « 7 <i)], 

where a"j refers to the a-component of the Pauli matrix at site j of the block labeled by /. The exact treatment of 
this Hamiltonian leads to four distinct eigenvalues which are doubly degenerate. The ground, first, second and third 
excited state energies have the following expressions in terms of the coupling constants. 

i^o) = -^L=(im>+9iut> + im» , w> = -i^(iTU>+9im> + iuT», 

e = -[ J 2 + A - \J {J 2 - A - A 2 f + 8], 

iv>i)= J_ (iffl)+PinT) + utt» > K)= J_ 0TU)+pim) + iuT)), 

V 2 + p 2 V 2 + p 2 



ex = -[ J 2 + A + \J {J 2 - A- A 2 f + ; 



i^) = ^(iut)-iru)) , m = ^(i m) - 1 m», 

e 2 = 2J 2 + A 2 ), 



w = im) . w> = uu>, 

e 3 = ^(A 2 -2A), 



where q,p are 



q = -^[J 2 - A - A 2 + y/(J 2 - A - A 2 f + 8] , p = -^[J 2 - A - A 2 - yj (J 2 - A - A 2 f + 8]. 
I |) and I J.) are the eigenstates of a z . 

In Fig. we have presented the different regions where the specified state is the ground state of the block 
Hamiltonian. The border between these regions are specified as a projection to a fixed plane. The projection to 
A 2 = J 2 A (A = 5) plane is shown in Fig[4ja) and the projection to A 2 = plane is plotted in FigfJJb). 
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APPENDIX B: THE RENORMALIZED COUPLING CONSTANTS OF THE EFFECTIVE 

HAMILTONIAN IN REGION A 
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APPENDIX C: THE RENORMALIZED COUPLING CONSTANTS OF THE EFFECTIVE 

HAMILTONIAN IN REGION B 
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